Large-deviation analysis for counting statistics in mesoscopic transports 
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We present an efficient approach, based on a number-conditioned master equation, for large- 
deviation analysis in mesoscopic transports. Beyond the conventional full-counting-statistics study, 
the large-deviation approach encodes complete information of both the typical trajectories and the 
rare ones, in terms of revealing a continuous change of the dynamical phase in trajectory space. 
The approach is illustrated with two examples: (i) transport through a single quantum dot, where 
we reveal the inhomogeneous distribution of trajectories in general case and find a particular scale 
invariance point in trajectory statistics; and (ii) transport through a double dots, where we find a 
dynamical phase transition between two distinct phases induced by the Coulomb correlation and 
quantum interference. 
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I. INTRODUCTION 

It has been widely recognized that, beyond the usual 
average current, the current fluctuations in mesoscopic 
transport can provide as well very useful information for 
the nature of transport mechanisms pj. In particular, 
a fascinating theoretical approach, known as full count- 
ing statistics (FCS)[2|, can conveniently yield all the sta- 
tistical cumulants of the number of transferred charges 
through a variety of mesoscopic systems (3l-[l3j. Experi- 
mentally, the real-time counting statistics has been real- 
ized in transport through quantum dots [l4j , representing 
a crucial achievement of being able to count individual 
tunneling events in mesoscopic transport. 

For mesoscopic transport, in addition to the well- 
known Landauer-Biittiker scattering theory |15| and 
the nonequilibrium Green's function formalism [16J, a 
particle-number-resolved master equation approach was 
demonstrated very useful in studying quantum noise in 
transport [l7rl2^ |. The key quantity in this approach 
is the particle-number(rt)-dependent density matrix (i.e., 
p( n \t), the reduced state of the central system in the 
transport setup). The connection of p^ n \t) with the 
distribution function in counting statistics simply reads 
P(n,t) = Tr[pW(t)], where the trace is over the central 
system states. From this distribution function, all orders 
of cumulants of the transmission electrons can be cal- 
culated, most conveniently, by introducing the cumulant 
generating function (CGF) e'^^ = J2 n P{n,t)e inx . 
This is in fact a discrete Fourier transformation, and x 
corresponds to the so-called counting field. 

Instead of the above discrete Fourier transformation, 
we may consider a replacement: e %xn => e~ xn . That 
is, we introduce the duaZ-function of P(n,t): P(x,t) = 
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J2 n e~ xn P(n, t) = e - - 7 ^ 2 ' 4 ). The real nature of the trans- 
forming factor e~ xn , in contrast to the complex one, e lxn , 
makes the resultant P{x, t) in some sense resemble the 
partition function in statistical mechanics. This function, 
in broader contexts, is named large- deviation (LD) func- 
tion, which, remarkably, allows to access the rare fluctua- 
tions or extreme events in random an d dy namical systems 
when performing the LD analysis [23l |. The partition 
function, in equilibrium statistical mechanics |24j . mea- 
sures the number of microscopic configurations accessible 
to the system under given conditions. The LD function, 
analogously, is a measure of the number of trajectories 
accessible for the dynamical aspects of a nonequilibrium 
system (23[. Relating to mesoscopic transports, if we are 
interested in the dynamical aspects of the transport elec- 
trons in time domain, again, analogous to the partition 
function, the LD function is a measure of the number of 
trajectories accessible to the electron "counter". Using 
the terminology in LD analysis, the trajectories in meso- 
scopic transports are categorized by the dynamical order 
parameter "n" or its conjugate field "x" . Trajectories are 
given an exponential weight, e~ xn , playing a role similar 
to the Boltzmann factor in statistical mechanics. 

Since the LD theory describes the trajectory space 
from multiple angles through the effect of the conjugate 
fields, say, analyzing the "x" dependence instead of lim- 
iting x ~> as in the FCS study, the LD approach is 
expected to be more instrumental in studying the statis- 
tics of quantum systems. In a most recent work [25| . 
the LD analysis was applied to a number of well-known 
simple quantum optical systems. Unexpectedly, striking 
dynamical phenomena such as scale invariance of trajec- 
tories as well as dynamical phase transition were uncov- 
ered for the statistics of photon emissions. These findings 
indicate that more surprising and complex dynamical be- 
haviors are seemingly waiting to be discovered in various 
scenarios. In the present work, we formulate a convenient 
LD method for mesoscopic transports and apply it to a 
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couple of examples. 

The paper is organized as follows. In Sec. II we 
begin with a brief outline of the particle-number (re- 
conditioned master equation together with a summary 
of the conventional FCS formalism based on it, before 
formulating the LD approach to quantum transport. As 
a first illustrative example, we present in Sec. Ill an LD 
analysis for the simplest system of transport through a 
single quantum dot where we reveal an inhomogeneous 
distribution of trajectories in general and find a partic- 
ular scale invariance point in trajectory statistics. We 
then carry out in Sec. IV the second example of transport 
through a parallel double dots where we find a dynamic 
phase transition between two distinct phases induced by 
the Coulomb correlation and quantum interference. Fi- 
nally, in Sec. V we summarize the work. 



t a)Jl the tunneling amplitudes under an approximation of 
energy independence. 

Remarks. — (i) The above master equation is up to 
the second-order cummulant expansion for the tunnel- 
ing Hamiltonian (26j . While this approximation applies 
well to most dissipative systems (e.g., in quantum op- 
tics), in quantum transport it actually limits the ap- 
plicability to the so-called sequential tunneling regime. 
Generalization to higher order expansion of H' is possi- 
ble, but quite complicated [27|- (H) The connection of 
the particle-number-conditioned density matrix with the 
distribution function in counting statistics is straightfor- 
ward, i.e., P(n,t) = Tr[p( n '(t)], where the trace is over 
the central system states. We will see in next subsec- 
tion that all the cumulants of current can be calculated 
conveniently from this distribution function. 



II. LARGE-DEVIATION METHOD 
A. Transport Master Equation 

In general, a mesoscopic transport setup can be de- 
scribed by the following Hamiltonian 

H = H s (al,a li )+ ^ J2e ak bl k b ak 

a=L,R k 

+ 2 I>«^a|> Qfe +H.c.). (1) 

a—L,R fik 

Here, Hs is the Hamiltonian of the central system in be- 
tween the electrodes, with aj^ (a M ) the creation (annihila- 
tion) operator of the (single particle) state u fx" . The sec- 
ond and third terms describe, respectively, the two (left 
and right) electrodes and the tunneling between them 
and the central system. 

Introducing F^ = Y, a k t afl kb a k = /l m + fa?, we re- 
express the tunneling Hamiltonian as H' = Y^u,( a n^ + 
H.c.) . Regarding the tunneling Hamiltonian as a pertur- 
bation, and properly clarifying the subspace of electrode 
states in association with the electron number transmit- 
ted through the system, say, the electron number ar- 
rived to the right electrode, we can obtain a number- 
conditioned master equation (l9| . which formally reads 

p [n) = - iC P {n) - ^ K jP {n+j) . (2) 

J=0,±1 

Here, p( n ' is the reduced density matrix of the system 
conditioned on the electron number "n" arrived to the 
right electrode. The Liouvillian C is the well-known com- 
mutator defined by the system Hamiltonian Hs ■ The su- 
peroperators IZj describe electron tunneling between the 
system and the electrodes, with an explicit form referred 
to Ref. . For the convenience of latter quotation, here 
we only mention an expression for the tunneling rates, 
based on the Fermi's golden rule, T a ^ = 2irtt a \t ai i\ 2 , 
where fl a are the density of states of the electrodes, and 



B. Counting Statistics 

With the knowledge of P(n, t) from the n-dependent 
density matrix p^ (t) , the current cumulant generating 
function (CGF) can be constructed through 

e -Hx,t) =J2P{n,t)e inx , (3) 

n 

where y is the so-called counting field. Based on the 
CGF, jF(x,t), the fcth cumulant can be readily carried 

out via Ck = —{-~idx) k ^ : {Xi't)\x=o- As a result, one 
can easily check: the first two cumulants, C± — n and 
C*2 = n 2 — n 2 , give rise to the mean value and the 
variance of the transmitted electrons, while the third 
cumulant (skewness), C3 = (n — h) 3 , characterizes the 
asymmetry of the distribution function. Here, (• ■ • ) = 
(• • • )P(n, t). Moreover, one can relate the cumulants 
to measurable transport quantities, e.g., the average cur- 
rent by I = eCi/t, and the zero-frequency shot noise by 
S = le^Cijt- Also, the important Fano factor is sim- 
ply given by F = C2/C1, which characterizes the extent 
of current fluctuations: F > 1 indicates a super-Poisson 
fluctuating behavior, while F < 1 a sub-Poisson process. 



C. Large-Deviation Formalism 

Instead of the discrete Fourier transformation, Eq. (|3|) , 
which is employed in the conventional counting statistics, 
let us now consider a replacement: e* x " => e~ xn . That 
is, we introduce the following dual- function of P(n,i): 

P(x, t)=J2 ^ xn P{n, t) = e"^'*). (4) 

n 

The real nature of the transforming factor e~ xn , in con- 
trast to the complex one, e txn , makes the resultant 
P(x, t) in some sense resemble the partition function in 
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statistical mechanics. Using analogous language in sta- 
tistical mechanics, in Eq. the trajectories are catego- 
rized by a dynamical order parameter "to" or its conju- 
gate field "x" . This is realized by the exponential weight 
similar to the Boltzmann factor, with the dynamical or- 
der parameter playing the role of energy or magnetiza- 
tion and the conjugate field the role of temperature or 
magnetic field. 

The function P(x,t), in broader contexts, coincides 
with the large-deviation (LD) function which plays an 
essential role in large-deviation analysis. In statistical 
mechanics, the partition function measures the number 
of microscopic configurations accessible to the system un- 
der given conditions. For the mesoscopic transport under 
consideration, if we are interested in the dynamical as- 
pects of the transport electrons, the above insight can 
lead to a LD analysis in time domain. That is, the LD 
function is a measure of the number of trajectories ac- 
cessible to the "counter", which favorably characterizes 
the trajectory space from multiple angles through the 
effect of the conjugate field. In particular, it allows to 
inspect the rare fluctuations or extreme events by tuning 
the conjugate field "x" . 

We emphasize that, if one performs only the conven- 
tional FCS analysis, using either a complex transform 
factor e %xn or a real one e~™ would make no difference 
since at the end the limit xi x ) ~ * wn l be involved. 
However, for the LD study, we must use the real factor 
e~ xn , which has a role of categorizing (selecting) the tra- 
jectories. It is right this type of selection which enables 
us to perform statistical analysis for the fluctuations of 
sub-ensembles of trajectories. For instance, x > implies 
mainly selecting the inactive trajectories (with small n), 
while x < prefers the active trajectories (with large 
n). In particular, by varying x, the x-dependent statis- 
tics can reveal interesting dynamical behaviors in time- 
domain. In other words, based on the distribution func- 
tion P(n, t) which contains complete information of all 
the trajectories, the LD approach, beyond the conven- 
tional FCS, captures more information from P(n,t) via 
the x-dependent cumulants. 

Let us now turn to the technical aspect of the LD ap- 
proach, based on Eq. ([2]). Similar to transforming P(n, t) 
to P(x,t), we introduce p(x,t) = J2 n e ~ xn P^ n \' t )- Then, 
from Eq. we obtain 

p(x, t) = [-iC -Ho- e x Ki - e- x TZ-i] p(x, t). (5) 

This equation allows us to carry out the LD function 
P(x,t), via P(x,t) = Tr[p(x,t)], where the trace is over 
the central system states. Accordingly, we obtain the 
generating function T(x,t) = — hiP(x,t), for arbitrary 
counting time t. Further, we can prove: 

= d x F{x,t) = -J—y2ne~ xn P(n,t) = (n) x , 

n 

(6a) 

F 2 {x,t)=d 2 x F{x,t) = -({n-n x )%, (6b) 



and more generally, J-k(x, t) = d^J r (x, t) = ( — )( fc +i) ((to— 
n x ) k )x- Here, for brevity, we used also the notation n x 
for (n) x . Using these cumulants, we can define a finitc- 
counting-time average current I(x,t) = eJ r i(x,t)/t, and 
the shot noise S(x,t) = 2e 2 J 7 2{x,t)/t. Also, the gener- 
alized Fano factor, F(x,t) = Tiix, t)/F\(x, t), will be of 
interest to characterize the fluctuation properties. 

We notice that, in order to obtain !Fk(x,t), we need 
only to carry out the various /c-th order derivatives of 
P(x,t), Pk(x,t) — d x P(x,t). This observation allows 
us an efficient method to compute the x-dependent cu- 
mulants for finite counting time. That is, performing 
the derivatives d k on Eq. ([5]) and defining pf.(x,t) = 
d x p(x,t), we get a set of coupled equations for pk(x,t), 
whose solution then gives Pk(x,t) = Tr\fik(x, t)]. 

In long counting time limit, it can be proved that 
T(x,t) ~ tX{x). In the remaining parts of this work, 
we will also call X(x) as an LD function, with no more 
distinction from P(x,t). Desirably, we notice that the 
asymptotic form, P(x,t) ~ e ~ tx ( x ) ; allows for a simpler 
way to get the LD function A(a;). That is, one can iden- 
tify it to the smallest eigenvalue of the counting matrix, 
given by the r.h.s of Eq. JS]). 

In this context, we would like to mention that the limit 
of x — > simply gives Fk{x, t)|x->o = Cfe, recovering the 
usual FCS cumulants. Again, we remark that the LD 
function around x = encodes information of the typical 
trajectories, while away from x = 0, on the other hand, 
it encodes information about the rare trajectories. For 
nonzero x, the possible anomalous behavior of X(x) with 
x, as we will see in the following examples, is usually 
a signature for nontrivial set of dynamical trajectories. 
For instance, singularities of X(x) correspond to certain 
dynamical (or space-time) phase transition, i.e., a rapid 
crossover between two distinct dynamical phases. 

III. TRANSPORT THROUGH A 
SINGLE-LEVEL QUANTUM DOT 

As a preliminary example, let us consider the transport 
through a single- level quantum dot [l4j|. Also, we assume 
spinless electron for simplicity but not losing physics in 
the sequential tunneling regime (large bias limit). In this 
case the system Hamiltonian in Eq. (|T|) has a simple form 
Hs = Eoa^a. which implies only two states involved in 
the transport process, i.e., the empty state |0) and the 
occupied one |1). Using this state basis, the n-dependent 
rate equation reads [19( 
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where R are the well known tunneling rates given by 
Ferm's golden-rule. Now, transforming Eq. (J7J in terms 
of p(x,t) = P (t)e~ xn and transposing p(x,t) into 
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a column-vector form [poo, Pu, Poi] T , we formally obtain 

p(x,t) = -M(x)p(x,t), (8) 

where Ai{x) is a 3 x 3 matrix, explicitly, with the follow- 
ing form 

I T L -e- x T R \ 

Mix) = -r L r R o . (9) 



For the purpose of an analytic analysis, we carry out 
the eigenvalues of Ai(x). They are, respectively, Ai = 
(r L + T r)/2, and A 2;3 = (r L + 1^ ± X)/2, where 
X = y/(T L -T R ) 2 +4e-*T L T R . Since X > 0, we 
can then identify that A3 is always the smallest eigen- 
value for arbitrary "x" , implying that the LD function 
X(x) = A3 (a). 

Knowing A(x), it is straightforward to calculate the x- 
dependent current and Fano factor. In Fig. 1 we plot 
A(x), I(x) and F(x) for the case of Tl = T R = T. 
We find that the x-dependent current decreases with x. 
This is because, through modifying the distribution func- 
tion P(n) to P(n)e~ nx , a larger x defines a less active 
subcnscmble of trajectories, which then gives a smaller 
subensemble average current. However, very interest- 
ingly, for this particular (equal coupling) configuration, 
we find in Fig. 1 an x- independent Fano factor, which 
indicates that all subensembles of trajectories, no matter 
how active or inactive, have the same fluctuation proper- 
ties of the typical trajectories. This remarkable scale in- 
variance of dynamical trajectories is very unusual, hold- 
ing only for T L = T R . 

Actually, we can carry out the analytic result of Fano 
factor for arbitrary coupling symmetry, which reads 
F{x) = l-\{\ + Re x )-\^eveR^{T L -T R ) 2 /AT L T R . 
Then, we see that, with the increase of x, the Fano fac- 
tor would increase and asymptotically reach the Pois- 
sonian limit of unity. This means that, differing from 
= T R , in the asymmetric setup the fluctuations 
will become stronger when we sample the trajectories 
from active to inactive ones. Finally, we notice that for 
highly asymmetric coupling, e.g., a = Tl/T r >> 1, 
F(x) ~ 1 - |(1+ f e*)- 1 ~ 1. This corresponds to a 
single barrier tunneling limit, which gives a Poissonian 
statistics for the tunneling events. 



IV. TRANSPORT THROUGH DOUBLE DOTS 
IN PARALLEL 

A. Model in a Transformed Representation 

In this section we consider a bit more complex setup 
which consists of a double (quantum) dots (DD) con- 
nected in parallel to two leads. Moreover, the inter- 
ference loop is pierced by a magnetic flux ($), similar 
to the well known Aharonov-Bohm interferometer. It 
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FIG. 1: (color online) The LD function X(x), s-dependent 
current I(x) and Fano factor F(x) for transport through a 
single-level quantum dot. In the symmetric case of Tl = 
F_r as assumed here, the x- independent Fano factor indicates 
a remarkable scale-invariance fluctuations of the dynamical 
trajectories. 



was shown recently that this system has some interest- 
ing properties such as a non-analytic current switching 
behavior, anomalous phase shift, and giant fluctuations 
of current — all of them induced by the interplay of inter- 
dot Coulomb correlation and quantum interference [28j . 

Following Ref. [HI, we assume that each dot has only 
one level, £-1(2)7 involved in the transport. In large bias 
limit, for simplicity, we also neglect the spin degrees of 
freedom, whose effect, under strong Coulomb blockade, 
can be easily restored by doub ling the tunneling rates 
of each dot with the left lead [28|. The entire system 
Hamiltonian reads 

H = H B + H T + E n a l a n + Ua\aia\a 2 . (10) 

/i=l,2 

a\ 2 are the creation operators for the DD, and the last 
term describes the interdot repulsion. The first term, 
h b = Y.k[ e LkV Lk b Lk + e R kb f Rk b Rk ], is for the leads and 
Ht for their couplings to the dots, 

H T = J2 (^Lalb Lk + t^b^a^ + H.c. . (11) 

Here, we assume that the couplings of the dots to the 
leads, i M L(fl)j arc independent of energy. In the absence 
of a magnetic field one can always choose the gauge in 
such a way that all couplings are real. In the presence of 
a magnetic flux $, however, the tunneling amplitudes be- 
tween the dots and the leads are in general complex. We 
write ^L(fl) = ^i(fl)e l< ^ L(H) , where is the cou- 

pling without the magnetic field. The phases are con- 
strained to satisfy 4>il + 4>iR ~ ^2L — 4>2 R = 4>, where 
<j> = 27r<f>/(f>o. In this work, we only consider the case with 
strong inter-dot Coulomb blockadc(£a, 2 + U ^> ^L)(two 
dots cannot be both occupied simultaneously). 

We found, in Rcfs. [2^|, that a transformed basis can 
benefit the analysis a lot. For simplicity, we consider 
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Ei = E<x- In this case, the SU(2) basis transformation 
reads 



a 2 



J_ I tlR t 2 R 

in 



L 2R b lR 



"1 

a-2 



(12) 



where tR = (t^ R + ^r) 1 ' 2 - Under this transforma- 
tion, the DD Hamiltonian, i^DD = X)u=i 2 E^a^a^ + 

U a\a\a\a2, is invariant, while the tunneling Hamiltonian 
is transformed to 



E 

jj. . k 



(t^a-lbkL + t^Rb^a^ + H.c. , (13) 



where, in the transformed basis, the tunneling ampli- 
tudes read 



tlR 
tlL 
<2Z 



tR, t 2 R — 0, 

e ^+t>™)(t lL t lR e^ + i2 L i 2R )/i R , 
_^ 2L - M ( iiL t 2R £ it _ t 2L t 1R )/tR. (14) 



matrix: 



M(x) 



( 2r L 


-2T R e 


-X Q 





°\ 


-(i + cos^)r L 


2T R 











-(i-cos0)r L 














—i sin^ 








Tr 





\ i sine/) Tl 











TrJ 



Its eigenvalues read, respectively, 

(r L + r R ) ± 



Ai = 0, A 



(16) 



and A 



2(3) 



4(5) 



X, where X 
^(r L -TR) 2 + 2e-xr L rR{l + cos<f>). Since A 2 , 3 ,4 are 
positive, we can conclude that the smallest eigenvalue 
is the smaller one among Ai and A5. The analysis in Sec. 
IV (D) will show that, by varying the conjugate field x, 
A5 has a crossover from negative to positive, implying 
that Ai = will suddenly become the smallest eigen- 
value. This behavior actually indicates a dynamic phase 
transition. Before proceeding to such an analysis, for the 
sake of completeness, in the following we first outline a 
scheme to account for the dephasing effect between the 
dots. 



The advantage of using the transformed basis is thus 
clear: the state aJlO) becomes decoupled with the right 
lead, and its coupling with the left lead is magnetic-flux 
tunable (in particular, can be switched off). 



B. Large-Deviation Function 

In the transformed basis, the strong Coulomb repul- 
sion among the inner dot and inter dots defines a reduced 
Hilbert space expanded by {|0), |1), |2)}, where |0) stands 
for the vacant DD state, and |1(2)) for the single occu- 
pation state a^ 2 )|0). Using this basis, the numbers- 
dependent master equation reads 



■ (« 
Poo 

.(n 
Pil 

.(n 
P22 

■(>' 
Pl2 

.{n 
P21 



= -2T LP ^+2TR P [r 1 \ 

(n) 



» 



=(i + cos0)r LjOo ^-2r^ii 
=(1 - cos0)r L ,4o , 

= -is\\\4>T L p y m ' - T R p 21 '. 



C. Dephasing Effect 

In real systems, dephasing may originate from sur- 
rounding environments, which cause random fluctuations 
of energy levels. For the present double-dot setup, a con- 
trollable dephasing mechanism can be introduced in ex- 
periments by performing a "which-path" detection using 
a nearby quantum point contact (QPC) [29l - l3lj |. The de- 
phasing effect of such a "which-path" detection can be 
accounted for by adding a Lindblad term [29| 



(») r.t 



1 



LL--L 



\L+pW - (17) 
The pure dephasing operator, 



(15) 



L 4>P 



to the r.h.s of Eq. (|I5 

La, 



\JT2(\l)r)(l\ — |2)ij(2|), describes the energy fluc- 
tuations of state |1)d and \2)d (in the original basis). 
Moreover, the dephasing rate can be carried out as |29| : 
Td = (VT - VT) 2 V d /27r, where V d is the bias voltage 
across the QPC, and T and T are the respective trans- 
mission probabilities, depending on which dot among the 
two is occupied. Note also that, since the dephasing does 
not affect electron's counting, the above dephasing term 
is not related to p^ n ~ x \ but only to p^ n \ 



Here, again, "n" is the electron number counted at the 
right lead. Also, for the sake of brevity, we assumed 
that iiL = i 2 L = i^, and t\R = I2R = tR. Then, the 
tunneling rates are defined, respectively, as Tl = 2^0^ 
and Tr = 2ir£lP R , with f2 the density of states of the 
leads. 

Similar to the single dot studied in the preced- 
ing section, making the LD transformation, p{x 1 i) = 
J2 n P^ n Ht)e~ xn , and transposing p(x,t) into a column- 
vector form [poo, P11, P22, P12, P2i] T , we can reexpress 
Eos. ([T5]) as p(x,t) = —M{x)p{x,t), with M(x) a 5 x 5 



D. Results and Discussions 

In the following, we restrict our numerical simulation 
to a symmetric setup, i.e., Yl = Tr = T. and use 
r (r" 1 ) as the units of energy (time). For the case 
of finite counting time (t c ), applying the scheme de- 
scribed in Sec. II (C), we identify the LD function from 
X(x,t c ) = F(x,t c )/t c and simply denote it as X(x). Sim- 
ilar simplified notation will be used also in the follow- 
ing for the a:-dependent ("x- ensemble") current I{x) and 
Fano factor F(x). 
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In Fig. 2 we display the LD function X(x) and the asso- 
ciated x-ensemble current I(x). We consider a relatively 
long counting time t c = lSOr -1 , and assume the dcphas- 
ing rate = 0.01T. Drastically differing from the above 
result of single dot, here, the LD function (solid lines in 
Fig. 2) reveals a phase transition behavior, when crossing 
through a ^-dependent critical point. 

To understand this behavior, we refer to the analytic 
solution obtained in Sec. IV (B), under the ideal case 
of Td — 0. Since the eigenvalues X2,3,i arc positive, the 
smallest eigenvalue, which has dominant contribution to 
the LD function in long counting-time limit, should be 
either Ai or A5, depending on the sign of A5. We can 
easily check that, as a function of the "dynamical field" 
x, A5 can be either positive or negative, with a crossover 
at x = 2 In I cos || = x c (<p). That is, for x > x c , A5 is 
positive, while for x < x c it is negative. This indicates 
that the LD function X(x) — for x > x c , and X(x) = 
X$(x) for x < x c . 




FIG. 2: (color online) The LD function X(x) (solid lines) and 
the i-dependent current I(x) (dashed lines) for a double-dot 
Aharonov-Bohm interferometer in strong Coulomb blockade 
regime. We assumed a symmetric setup with Fl = Fr — F 
and Ei = E2, a relatively long counting time t c = 150F - , 
and a dephasing rate Fa — O.Oir. 



As a further complementary analysis, let us evaluate 
the distribution function P(n, t) in long time limit. Ow- 
ing to the assumed Coulomb blockade effect, we know 
that once the (transformed) upper dot is occupied, i.e., 
the DD is in the state ^l^), the current would be 
completely blocked. Denoting the tunneling probabil- 
ity from the left electrode to the upper (lower) dot 
as P2(i)i the probability with "n" electrons transmitted 
through the DD (actually through the lower dot) sim- 
ply reads P{n) = (pi) n p2, before the transport chan- 
nel is blocked. This distribution function defines the 
entire ensemble of trajectories, and supports as well a 
straightforward "subcnsemble" LD analysis. Noting that 
Pi(2) =_ f i(2)l(0)/[? il{4>) + r 2 i (0)], for the_symmetric 
setup t\L = hh = t L and ti R = t 2 R = t R , we ob- 
tain a simple result: p\ = cos 2 (^), and p 2 = sin 2 (^). 
Now, let us consider the LD transformation, P{x) = 
p 2 X)^Lo e ~ xn (pi) n ■ This infinite series has a convergent 
sum, P(x) = sin 2 (^)[l — e~ x cos 2 (^)] _1 , under the con- 



dition e x cos 2 (|) < 1. This convergence condition gives 
x > 21n| cos(f )|, implying that the LD function (A(x)) is 
zero in this region, being fully consistent with the result 
based on the eigenvalue analysis. 

In Fig. 2, we also plot the £-dependent current I{x), by 
the solid lines. This re-sampled current over a subcnsem- 
ble of trajectories shows a switching behavior around the 
particular point x c (<ft). In the ideal case of E\ = E%, 
Td = and t c — > 00, the switching behavior of I(x) 
deforms to a discontinuous jump, owing to the singu- 
larity of X(x) at x c {4>)- This feature clearly indicates a 
first-order dynamical phase transition, revealing a rapid 
crossover between two distinct dynamical phases across 
x c : an active phase given by x < x c with trajectories 
having large n, and an inactive phase by x > x c with 
trajectories having small n. At the critical point x c (<j)), 
the two dynamical phases coexist. 
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FIG. 3: (color online) Effects of the counting time t c and 
the dot-level offset A. In (A), for a given A = O.ir and 
4> = 0.25-7T, we display the results of t c = 150r _1 (blue dotted 
lines) and 50r -1 (red dashed lines) against that of the long t c 
limit (black solid lines). In (B), in the limits of long counting 
time and <f> = 0, we compare the results of A = O.ir (black 
solid lines), 0.5r (blue dotted lines) and F (red dashed lines). 
In both (A) and (B), we assumed Fl = T_r = F and Fd ~ 0. 
The latter assumption is mainly for the purpose of showing 
the pure effects of the counting time and the level offset. 



In Fig. 3 (A) and (B) we display, respectively, the ef- 
fects of the counting time (t c ) and the dot-level offset 
(A). In Fig. 3 (A), as a comparison to the result of long 
t c limit (black solid lines), we plot the LD function, the 
current and the Fano factor for two finite counting times, 
by the blue dotted and red dashed lines. We sec that, as 
expected, a shorter counting time will result in a milder 
crossover behavior. As a similar comparison, in Fig. 3(B) 
we show the effect of the dot-level offset. Since in the 
transformed basis the two dots are coupled with a tun- 
neling amplitude proportional to A, we then understand 
that the dynamical phase transition will be less drastic 
as we increase A. 
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V. SUMMARY 

Starting with a number-conditioned master equation, 
we formulated an efficient large-deviation approach for 
mesoscopic transports. Formally, it appears that the 
large-deviation approach is similar to the FCS. However, 
a large-deviation analysis is capable of capturing more 
complete information: not only for the typical trajecto- 
ries, but also for the rare ones. Through an inspection 
on the continuous change of subensemble statistics of 
trajectories, richer dynamical behaviors can be revealed. 
Methodologically, in addition to the eigenvalues analysis 
in long time limit, we also presented a convenient scheme 
of calculating the large-deviation function for arbitrary 
counting times. 

As illustrative applications, we considered two exam- 
ples: one is the transport through a single dot, another 
is through a double dots. Even for the simple system 
of single dot, we found that the symmetry of the tunnel 
coupling to the leads would strongly affect the statistics 
of the trajectories, which is indeed outside the scope of 
the usual FCS analysis. Very interestingly, similar to 



the "unexpected" finding uncovered in the quantum op- 
tical system of driving two-level atoms J2f| , here a sym- 
metric tunnel-coupling in our single-dot transport sys- 
tem can as well result in a ("time" or "tunneling-rate" ) 
scale-invariance behavior in the statistics of trajectories. 
For the double-dot setup, our large-deviation analysis re- 
veals a more profound behavior, say, the dynamical phase 
transition, which is, essentially, induced by the Coulomb 
correlation and quantum interference, and were analyzed 
in detail by combining numerical simulations with some 
analytic solutions. 

The large-deviation approach, owing to its ability to 
analyze the subensemble statistics of trajectories, is an- 
ticipated to apply to various scenarios, e.g., from elec- 
tronic nano-devices to femtochemistry. In particular, its 
systematic applications to important mesoscopic trans- 
port systems are extremely promising, which is seemingly 
a new area waiting for exploitation. 
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